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ABSTRACT 

Hills (1988) predicted that runaway stars could be accelerated to velocities larger than 
1000 km s" 1 by dynamical encounters with the supermassive black hole (SMBH) in the 
Galactic center. The recently discovered hypervelocity star SDSS J090745. 0+024507 
(hereafter HVS) is escaping the Galaxy at high speed and could be the first object in 
this class. With the measured radial velocity and the estimated distance to the HVS, 
we trace back its trajectory in the Galactic potential. Assuming it was ejected from the 
center, we find that a ~ 2 mas yr _1 proper motion is necessary for the star to have come 
within a few parsecs of the SMBH. We perform three-body scattering experiments to 
constrain the progenitor encounter which accelerated the HVS. As proposed by Yu 
& Tremaine (2003), we consider the tidal disruption of binary systems by the SMBH 
and the encounter between a star and a binary black hole, as well as an alternative 
scenario involving intermediate mass black holes. We find that the tidal disruption of 
a stellar binary ejects stars with a larger velocity compared to the encounter between 
a single star and a binary black hole, but has a somewhat smaller ejection rate due to 
the greater availability of single stars. 
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1 INTRODUCTION 



The recent discovery bv lBrown et al.l J2005) of the hyperve- 
locity star SDSS J090745.0+02450 7 (he r eafter the HVS) has 
lent credence to the prediction of iHillsl il988l) that dynam- 
ical encounters with the black hole in the Galactic centre 
can eject stars with velocities up to several thousand kilo- 
meters per second. At a distance of 40-70 kpc, the HVS has 
a heliocentric radial velocity of 853 ± 1 2kms -1 in a direc - 
tion of 174° from the Galactic centre ferown et alJ 12005). 
Once corrected for solar motion and galactic rotation, this 
translates into a velocity of about 730kms _1 relative to the 
Local Standard of Rest. This velocity, which represents the 
lower limit of the total space velocity, is significantly higher 
than th at of any oth er runaway or high-velocity star in the 
Galaxy dStonelll99ih . 

We investigate the origin of the extreme velocity of the 
HVS by means of kinematic analyses, binary evolution cal- 
culations and numerical scattering experiments. There are 
three possible explanations for the velocity of the HVS: (i) 
ejection upon supernova explosion of the companion in a bi- 
nary system; (ii) dynamical ejection after a close encounter 
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with main sequence stars or stellar mass compact objects; 
(iii) dynamical ejection from the Galactic centre as a result 
of an encounter with the supermassive black hole. 

The paper is organised as follows: in §2 we analyse the 
possibility of an origin in the Galactic disk and an ejec- 
tion by supernova explosion, in §3 we back-trace the orbit 
of the star in the Galactic potential and constrain the value 
of its proper motion and in §4 we describe numerical scat- 
tering experiments of 3-body encounters involving main se- 
quence (MS) stars, the supermassive black hole (SMBH) and 
intermediate-mass black holes (IMBHs). 



2 THE ORIGIN OF THE HIGH VELOCITY OF 
THE HVS 

To address the possibility that the HVS is a high-velocity 
runaway resulting from the disintegration of a binary sys- 
tem, we turned to population synthesis. Making use of th e 
SeBa stellar evolution package iPorteeies Zwart et al.l200lf) . 
we generated two sets of 10 6 binaries each. In the first set, 
primary masses ranged from 100 Mq down to the hydro- 
gen burning limit. In the second set, the minimum primary 
mass was increased to 8Mq, in order to get a larger sam- 
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pie of events for statistical purposes. The mass ratio i n bot h 
cases was chosen from the distribution of iHoeeveenl il992f) . 
whilst the initial orbi t al sep aration was drawn from that of 
iDuauennov fc Mavorl (ll99lT) and truncated at 10 5 R© . The 
distr ibution of natal k ick speeds for neutron stars was taken 
from IPaczvnskil l|l990l l. with a dispersion of 300kms _1 . We 
assumed that the escape speed was the orbital speed of the 
secondary immediately before the binary disintegrated. As 
well, we only considered secondary stars with a mass below 
10 Mq, given the known constraints on the mass of the HVS. 
From this synthesis experiment, we obtained maximum es- 
cape speeds of ~ 70kms _1 , an order of magnitude below 
the space velocity of the HVS. Based upon this, we reject 
the high-velocity runaway hypothesis. The possibility that 
the system is a binary and that its high speed is the result 
of an asymmetric supernova kick can be discounted because, 
if we assume that the visible companion is on the main se- 
quence, a kick magnitude in excess of 2000 km s" 1 would 
be required. While this can not be ruled out physically, the 
probability of the system remaining bound in such an event 
is negligible. 

The only alternative explanation for the high veloc- 
ity of the HVS is an ejection during a dynamical en- 
counter. Stellar encounters involving MS stars, neutron 
stars or stellar mass black holes eject stars with maxi- 
mum velocitie s on the order of the orbital velocity in bi- 
nary systems ij Gruahmdrjj^^orto Eggletonl l2004t 
ISieurdsson fc Phinnevl Il993ft . iBrown et all J2005I) assume 
therefore that the HVS was ejected from the Galactic cen- 
tre by a strong encou nter with th e SMBH. This possibility 
was first proposed by [Hilla (Il98c1) . who considered encoun- 
ter s of binaries wit h a SM BH, and then further developed 
bv lYu fc Tremainel <l2003l) . who found that encounters be- 
tween binaries and the SMBH or between single stars and 
a binary black hole represent the most efficient channel to 
eject hypervelocity stars. 
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Figure 1. Trajectory of the HVS integrated backward in time 
in the Galactic potential without any proper motion component 
(dashed line) and with a proper motion of about 1.8masyr — 1 
(solid line). The integration is stopped upon passage through the 
disc. The full dot represents the present position of the star at 
a distance of 55 kpc while the grey dots represent a schematic 
model of the Galaxy. 



Table 1. Kinematic parameters of the HVS which minimize the 
distance from the Galactic centre. For each distance value, we 
report the average proper motion in right ascension fi a and dec- 
lination fig such that d m i n <10pc, the corresponding heliocentric 
velocity Vsun an d the velocity V e j at a distance d m ; n from the 
centre corrected for galactic rotation. 



d 


A*ct 




Vsun 




(kpc) 


(mas/yr) 


(mas/yr) 


(km/s) 


(km/s) 


39 


-1.4916 ±0.0007 


2.1679 ±0.0009 


982 ±1 


1246 ± 144 


55 


-0.9769 ± 0.0004 


1.5379 ±0.0003 


976 ±1 


1266 ± 66 


71 


-0.7214 ±0.0002 


1.1917 ±0.0003 


973 ±1 


1249 ± 142 



3 THE TRAJECTORY OF THE HVS IN THE 
GALAXY 

In the previous section we argue that the HVS must have 
been ejected from the Galactic centre by a dynamical en- 
counter with the SMBH. This implies that the HVS orig- 
inated in the Milky Way's central region. Despite the in- 
completeness of available kinematic data, we trace back the 
orbit of the HVS in the Galactic potential. The distance of 
the object is still uncertain , as it depends on its spectral 
type and evolutionary state. IBrown et alJ <l2005l) estimate a 
distance of 71 kpc for a B9.2 MS star and 39 kpc for a blue 
horizontal branch star, with an average value d = 55 kpc. 
We consider all three values in our analysis. We first assume 
that the star has no proper motion and trace its trajectory 
backward in time until it crosses the disc, using Paczyn- 
ski's model (1990) for the potential of the Galaxy. The or- 
bit, shown in Fig. (dashed line) for the case d = 55 kpc, 
doesn't pass through the Galactic centre region. We then 
randomly generate the two proper motion components in 
the range l-3masyr _1 and derive the three-dimensional po- 
sitions and velocities in the galactocentric reference frame 
J Johnson fc Soderbloml Il987l) . For each distance value we 
generate about 5000 sets of initial conditions, integrate the 



trajectories and determine the minimum distance d m in to 
the Galactic centre. We find that there is at least one com- 
bination of proper motion components for each distance such 
that d m in < 5pc. In Tab. we report, for each distance d, 
the average values /j, a and \i& of the proper motion compo- 
nents such that d m in <10pc, the corresponding heliocentric 
velocity Vsun and the velocity V c j at a distance d m in from 
the centre, corrected for galactic rotation. The ejection ve- 
locity V e j is about 1250kms _1 , independent of the assumed 
distance. In the case d = 55 kpc, the mean values fi a and 
/ia given in Tab. Q result in the orbit shown in Fig. Q (solid 
line) . 

We can use this analysis to calculate the minimum ve- 
locity with which the HVS should have been ejected from 
the disk in the case of the binary supernova scenario. As- 
suming a random proper motion in the range l-3masyr _1 , 
the minimum velocity at disk crossing (corrected for galactic 
rotation) is about 500kms _1 , much larger than the average 
recoil velocity from a supernova explosion (see § 2). This re- 
sult supports the scenario of an origin in the Galactic centre. 
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Figure 2. Examples of three-body encounters involving the SMBH. The positions are in units of the initial binary semi-major axis. 
(Left) An encounter between a binary of 3 Mq MS stars and the SMBH resulting in the the ejection of one star while the other star 
remains bound to the SMBH (exchange). (Middle) An encounter between a star and a black hole binary resulting in a preservation. 
(Right) An encounter between the SMBH and a binary containing a 3 Mg MS star and an IMBH resulting in the breakup of the binary 
and the ejection of the MS star (ionization). 



4 THREE-BODY SCATTERINGS WITH THE 
SUPERMASSIVE BLACK HOLE 

We now explore the hypothesis of a dynamical ejection 
from the Galactic centre by means of numerical simula- 
tions of three-body scatterings with the SMBH. In Fig.|5]wc 
show three examples of encounters involving MS stars, the 
SMBH and an IMBH. The experiments are carried out with 
the sigma3 package included in the STARLAB 1 software 
envir onment iMcMillan fc Hii3ll996t iPortegies Zwart et alJ 
l200ll) . For each simulation we specify: the masses of the 
three stars, the semi-major axis and eccentricity of the bi- 
nary and the relative velocity at infinity between the bi- 
nary's centre of mass and the single star. In order to classify 
possible collisions and mergers, physical radii are also spec- 
ified for the stars. Additional parameters like the orbital 
phase of the binary and its orientation relative to the in- 
coming star are rando mly drawn from uniform distributions 
<Hut fc Bahcalll985t . The initial ecc entricity is drawn from 
a thermal distribution llHeggielll975l) . The impact parame- 
ter b is randomized according to an equal probability distri- 
bution for b 2 in the range [0 — 6 max ]. The maximum value 
fcmax is determined automatically for each ex periment (see 
iGualandris. Portegies Zwart. Eggletonl J2004f) for a descrip- 
tion). Energy conservation is usually better than one part 
m 10 6 and, in case the error exceeds 10 5 , the encounter is 
rejected. The accuracy in the integrator is chosen in such a 
way that at most 5% of the encounters are rejected. 

In all the experiments, we consi der MS sta rs of mass 
m = 3 Mq (as indicated by iBrown et alJ <l2005l) for a B9 
star) and radius R = 2.4 R©, an IMBH of mass ibbh = 
3000 M w and a SMBH of mass Mm = 3.5 x 10 6 M 
jGhez et al J 120031 : ISchodel et al]l2003l) . The relative veloc- 
ity at infinity between the single star and the binary's cen- 
tre of mass is set equal to the dispersion velocity in the 
Galactic centre (~ 100 kms -1 ). The ejection velocities of 
escapers are taken at the distance at which the integrator 
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stops dMcMillan fc HutJll996l) . The effect of the black hole 
potential is negligible at this distance. 

4.1 Encounters between a binary of MS stars and 
the SMBH 

Close encounters between a binary and a very massive object 
can (i) break up the binaries and eject the two components 
with high speed (ionizations), (ii) eject one star and leave 
the second star bound to the SMBH (exchange) (see the left 
panel of Fig.|2Jl. We perform scattering experiments between 
the SMBH and binaries of MS stars. In all the runs, the 
binary stars have equal mass m and the semi-major axis is 
varied in the range 0.05 AU< a < 1 AU. 

In Fig. |3]we show the average ejection velocity Vcj of 
escapers as a function of the initial binary semi-major axis. 
Sufficiently high ejection velocities (V c j > 1250 kms -1 ) are 
obtained for a< 0.3 AU. The dotted l i ne rep resents the theo- 
retical prediction bv lYu fc Tremainel i2003l) with an ejection 
speed parameter = 130 kms -1 (see Eq. 20). 

In Fig. 3] we show the fraction of encounters resulting 
in ionization, exchange of the secondary star or merger for 
the range of a under consideration. The fraction of encoun- 
ters resulting in an exchange of the secondary star increases 
slowly from about 25% for a = 0.005 AU to about 30% for 
a — 1 AU. Since the binary components have equal masses, 
the probability of ejection for the two stars is equal in ex- 
change encounters. Ionizations are rare and only occur for 
a> 0.2 AU, w hich marks t he transition between hard and 
soft binaries iHeggid|l975l) . As a result, high ejection ve- 
locities are mostly produced in exchange encounters. The 
total fraction of scatterings whose outcome is a star escap- 
ing with a velocity larger than 1250 kms -1 decreases rapidly 
with increasing a, and therefore binaries with semi-major 
axis 0.05 AU <a< 0.3 AU are the most suitable progenitors 
of hypervelocity stars. If we consider that about 20% of the 
simulated encounters result in the ejection of a hypervelocity 
star over the suggeste d range of orbital separ ations, the ejec- 
tion rate provided bv lYu fc Tremainel i2003T) can be refined 
to ~ 3 x 10 -6 (77/O.I) yr" 1 , where 77 is the binary fraction. 



4 Gualandris, Portegies Zwart, Sipior 



10 4 



-T u 



□ □ □ 



1000 



n 



A A 



100 



(m,m) + SMBH 
SMBH = 3.5x10" M s 

m = 3 M„ 



10 



0.1 



a (AU) 



10" 



1000 



100 



1 1 

(SMBH, 


IMBH) + 


1 : 
m ! 


SMBH = 


= 3.5x10 s 


m s : 


■ IMBH = 


3000 M s 


- 


• 






• 






• 






• 






• 






i 






• 
















• 






i 



10 100 1000 

a (AU) 



Figure 3. Average recoil velocity of escapers as a function of 
the initial binary semi-major axis in the interaction of a stellar 
binary with the SMBH. The error bars indicate the 2a deviation 
from the mean. The squares indicate the velocity V m ax for which 
1% of the encounters have V c j > Vmax- The triangles indicate the 
velocity V m i n for which 1% of the encounters have V e j < V m i n . 
The horizontal line marks the 1250 km s — 1 ejection velocity of 
the H VS while t he do tted line gives the theoretical estimate by 
lYu fc Tremainj l2003t) . 
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Figure 4. Branching ratios as a function of the initial binary 
semi-major axis for encounters between a stellar binary and the 
SMBH. The empty symbols indicate the total fraction of ioniza- 
tions, exchanges and mergers while the full dots indicate the frac- 
tion of encounters which produce escaping stars with a velocity 
larger than 1250 km s — 1 . 

The fraction of physical collisions and mergers decreases 
steadily with increasing a from ~20% to ~2%, and mainly 
involve the binary components. Collisions with the SMBH 
occur in less than 1% of the cases. The merger products 
can remain bound to the SMBH or escape its gravitational 
potential. We perform additional scattering experiments of 



Figure 5. Average recoil velocity of a single MS star escaping 
from a black hole binary as a function of the binary's initial semi- 
major axis. The horizontal line marks the 1250 km s" 1 ejection 
velocity of the HVS while the dotte d line represents the analytical 
estimate bv lYu &: Tremaind i2003t) . 



encounters between binaries and the SMBH to study the 
properties of the merger products. We consider equal mass 
binaries with mass m = 1.5 Mq (in such a way that the mass 
of any possible merger is ~ 3Mq) and radius R — 1.4 R©. 
Escapers are ejected with velocities larger than 1000 km s -1 
if the initial semi-major axis is in the range 0.03 AU <a< 
0.05 AU. In this range, only about 6% of the encounters re- 
sult in a binary merger with escape of the collision product. 
We therefore consider this scenario inefficient for the pro- 
duction of hypervelocity stars. 



4.2 Encounters between a single star and a binary 
black hole 

We now consider the hypothesis that the SMBH is in a bi- 
nary with an IMBH and interacts with single stars (see the 
central panel of Fig. |5J . The semi-major axis of the black 
hole binary is taken in the range 2AU < a < 1000 AU. All 
the simulated encounters result in a preservation of the black 
hole binary, during which the single star gains kinetic energy 
and escapes. 

In Fig. |^| we show the average ejection velocity of 
the escaping star as a function of a. The maximum ve- 
locity obtained in these encounters is about lOOOkms -1 , 
barely sufficient to explain the velocity of the HVS. The 
comparison with the theoretical estimate (dotted line) by 
lYu fc Tremaind i2003t) reveals a discrepancy of about a fac- 
tor of 2 in the ejection velocities. The numerical results also 
flatten at large initial semi-major axes and deviate from the 
expected V c j oc a _1,/2 scaling. We have performed additional 
calculations with 1 Mq incoming stars but the average U c j 
appears to be insensitive to the mass of the single star. 
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Figure 6. Average recoil velocity of escapers as a function of 
the initial binary semi-major axis in the interaction between a 
normal star in a binary with an IMBH and the SMBH. The error 
bars indicate a 2a deviation from the mean. The squares and 
triangles are defined as in Fig. [3] The horizontal line marks the 
1250 km s — 1 ejection velocity of the H VS while the do tted line 
gives the theoretical estimate bv lYu fa Tremainel i2003l) . 

4.3 Encounters between a SMBH and a main 
sequence star in a binary with an IMBH 

It has been proposed dHansen fc Milosavlievic| I2003T) that 
the Galactic centre is populated by IMBHs, with masses 
in the range 100 - 10000 M Q . These IMBHs may form in 
young dense star clusters as a result of runaway collisions 
iPortegies Zwart et al.l I2004T) and sink toward the Galactic 
centre together with their parent cluster due to dynamical 
friction. While the cluster dissolves in the Galactic poten- 
tial, the IMBHs continue to spiral in, possibly with a stellar 
companion, until they eventually interact with the SMBH 
febisuzaki et alJl200ll) . If this is the case, IMBHs must play 
a role in the dynamical encounters taking place in the few in- 
ner parsecs of the Galaxy. We consider encounters between 
the SMBH and a binary consisting of a MS star and an 
IMBH. We select the binary semi-major axis in the range 
0.1 AU< a < 100 AU, with the lower limit set by the IMBH's 
tidal radius Rt = R (mBn/m) 1 ^ 3 = 24 R.0 for the adopted 
masses and radii. We focus on the encounters whose final 
outcome is the break-up of the binary with subsequent ejec- 
tion of star m to infinity, while the IMBH can either remain 
bound to the SMBH or escape (see the right panel of Fig. 

Figure HJ reports the velocity of the escaping star after 
the encounter as a function of the initial a. This type of en- 
counter can easily eject stars with velocities of thousands of 
kilometers per second . A theoretical estimate derived from 
lYu fc Tremainel fcOOSt (see Eq. 1-3) is shown with a dotted 
line. The numerical results agree well with the theoretical 
estimate. 

In Fig. |7Jwe show the fraction of encounters resulting in 
ionization, exchange of the secondary star or merger for the 
range of a under consideration. For a< 0.3 AU the binary is 
hard and the total cross-section is dominated by exchange 
encounters. In this case, the escaping star gains energy at the 



Figure 7. Branching ratios as a function of the initial binary 
semi-major axis in the interaction between a normal star in a 
binary with an IMBH and the SMBH. 

expense of the IMBH, which becomes bound to the SMBH. 
For larger semi-major axes, the binary is soft and tends to be 
ionized. The highest recoil velocities are generally obtained 
in exchange encounters. 

Although this scenario can eject stars with hyperveloc- 
ities over a wide range of orbital separations, only rarely ec- 
centric orbits result in prompt ionization. Prompt ionization 
only occurs for low angular momentum orbits, which com- 
prise a fraction of about 10 -5 if we assume an isotropic veloc- 
ity distribu tion. Adopting a form a tion r ate of ~ 10 _7 yr _1 
for IMBHs IPortegies Zwart et alJ i2005t) . the ejection rate 
can be as small as ?7l0 -11 yr -1 , where r] represents the frac- 
tion of IMBHs with a stellar companion. Low eccentricity 
orbits result in a slow inspiral of the binary due to dynam- 
ical friction, but, in this case, it is not clear whether the 
stellar density is sufficiently high to drag the binary to the 
tidal radius of the SMBH. If an encounter with the SMBH 
does take place, we expect the ejection velocity to be much 
lower than in the case of a prompt ionization. 



5 SUMMARY AND CONCLUSIONS 

We have investigated the origin of the extreme velocity of the 
HVS by means of kinematic analyses, binary evolution calcu- 
lations and numerical simulations of three-body encounters. 

By tracing the trajectory of the HVS in the Galactic 
potential (using the available measurements for the distance 
and the radial velocity) , we have shown that a proper motion 
of about ~ 2masyr _1 is required for the star to have come 
within a few parsecs from the SMBH. 

We confirm the prediction by iHillsl lll98ct) and 
lYu fc Tremainel i2003T) that dynamical encounters with the 
SMBH can eject hypervelocity stars to the Galactic halo; 
the HVS is likely the first discovered object of this kind. The 
most promising scenarios are dynamical encounters of MS 
stars with the SMBH, possibly in a binary with an IMBH. 
In particular, the encounter between a SMBH and a stellar 
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binary or between a single star and a binary black hole pro- 
vide enough kinetic energy to eject hypervelocity stars over 
a restricted range of initial semi-major axes. We have also 
investigated the more exotic encounter between a SMBH 
and an IMBH orbited by a stellar companion. Although this 
type of encounter can be very energetic and doesn't require 
any constraint on the binary's hardness, it has a very low 
probability. The possibility that the HVS is the product of 
a merger induced by the SMBH is intriguing but not very 
likely. 
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